Astron. Nachr. / AN , No. , 1 -[K)](2007) / DOI please setDOI! 



Combining astrometry with the light-time effect: 
The case of VW Cep, C Phe and HT Vir. * 

P. Zasche** and M. Wolf 

Astronomical Institute, Faculty of Mathematics and Physics, Charles University Prague, CZ-180 00 Praha 8, V 
Holesovickach 2, Czech Republic 

Received Today, accepted Today 
Published online later 



Key words binaries: eclipsing - stars: fundamental parameters - stars: individual (VW Cep, Zeta Phe, HT Vir) 

Three eclipsing binary systems with astrometric orbit have been studied. For a detailed analysis two circular-orbit binaries 
(VW Cep and HT Vir) and one binary with an eccentric orbit (C, Phe) have been chosen. Merging together astrometry and 
the analysis of the times of minima, one is able to describe the orbit of such a system completely. The O — C diagrams 
and the astrometric orbits of the third bodies were analysed simultaneously for these three systems by the least-squares 
method. The introduced algorithm is useful and powerful, but also time consuming, due to many parameters which one is 
trying to derive. The new orbits for the third bodies in these systems were found with periods 30, 221, and 261 yr, and 
eccentricities 0.63, 0.37, and 0.64 for VW Cep, £ Phe, and HT Vir, respectively. Also an independent approach to compute 
the distances to these systems was used. The use of this algorithm to VW Cep gave the distance d — (27.90 ± 0.29) pc, 
which is in excellent agreement with the previous Hipparcos result. 
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1 Introduction 

During the last century, many observations of close binary 
systems were collected, especially for eclipsing binaries. 
Most observations of these stars were made photometrically 
but in some cases also the spectroscopy was obtained. In a 
few cases, also astrometric observations were carried out 
and these systems, which were studied by several different 
techniques, are the most interesting, because one is able to 
find more relevant parameters of them. Especially, if it is 
possible to analyse the different measurements together (it 
means from different observational techniques in one least- 
squares fit) one is able to find the complete set of parame- 
ters describing the orbit of such a binary, which are not in 
contradiction. Such approach is potentially very powerful, 
especially in upcoming astrometric and photometric space 
missions. 

There are a few well-known eclipsing binaries with as- 
trometric measurements for which the light-time effect (he- 
reafter LITE) was considered, or expected. However, the 
eclipsing binary (hereafter EB) nature and the astrometric 
orbit were usually studied separately. 

For example, there are several published light-curve so- 
lutions for V505 Sgr with the third light included (see e.g. 
Lazaro et al. (2006), ibanoglu et al. (2000), etc), and an 
analysis of its apparent orbital period changes interpreted 
as the LITE due to the third body (e.g. Rovithis-Livaniou 
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et al. 1991). The only paper which compares the astrometry 
and a period analysis of O — C deviations from the constant 
orbital period was published by Mayer (1997). Despite ex- 
isting astrometric measurements, there were no attempts to 
combine these two methods together. The results from dif- 
ferent approaches were just compared to each other. 

Another systems are for instance QS Aql, 44i Boo, 
QZ Car, SZ Cam, GT Mus, or V2388 Oph. The coverage 
of the astrometric orbit is very poor for some of them. In 
the system QS Aql the LITE could be determined precisely, 
but the astrometric orbit is covered only very poorly, see e.g. 
Mayer (2004). The opposite case is V2388 Oph, where the 
astrometric-orbit parameters were computed very precisely, 
but there are only a few minimum-time measurements (see 
e.g. Yakut et al. 2004). For SZ Cam only a few usable as- 
trometric observations were obtained, but the LITE is well- 
defined and also the third light in the light-curve solution 
was detected (Lorenz et al. 1998). QZ Car is a more compli- 
cated, probably quadruple system consisting of an eclipsing 
and a non-eclipsing binary, but there are also only a few us- 
able astrometric measurements. Also GT Mus is a quadru- 
ple system, consisting of an eclipsing and RS CVn compo- 
nent. In many other cases, only measurements of the times 
of minima are available, without astrometry. For some oth- 
ers, astrometry without photometry, is only available. Other 
systems where the astrometric observations were obtained 
and the LITE is observable or expected are listed in Mayer 
(2004). 

The only paper on combining these two different ap- 
proaches into one joint solution is that by Ribas et al. (2002), 
where a similar method as described in this paper was ap- 
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plied to the system R CMa, but where only a small arc of the 
astrometric orbit was available. Besides the astrometry and 
LITE also the proper motion on the long orbit was analyzed. 
On the other hand one has to note, that in Ribas et al. (2002) 
the complete astrometric parameters (with proper motions, 
parallax, etc.) were used, while in this paper only a relative 
astrometry of the distant body relative to the eclipsing pair 
was analysed. 

2 Methods 
2.1 Astrometry 

The number of visual binaries with astrometric orbits has 
grown, but the whole orbit is covered with data points only 
in some cases. Thanks to precise interferometry the observ- 
able semimajor axes of astrometric binaries are still decreas- 
ing down to milli- and micro-arcseconds.On the other hand, 
one has to regret that no recent astrometric measurement 
of a wide pair of about 1" has been obtained for the sys- 
tems mentioned below. All the astrometric observations 
were adopted from 'The Washington Double Star Catalo- 
gue' WDS3 (Mason et al. 2001). 

From astrometric data (N measurements of the position 
angle Oi, separation pi, the uncertainties cr$, and <j p and the 
time of the observation ti) one is trying to find out the pa- 
rameters of the orbit for the distant component, defined by 7 
parameters: period P3, angular semimajor axis a, inclination 
i, eccentricity e, longitude of the periastron 073, the longi- 
tude of the ascending node il, and the time of the periastron 
passage T . One has to solve the inverse problem 

{{ti, 6i, Pi,cj0 ti ,a p ^)} l= i tN — > (a,p 3 , i, e, lu 3 , fl, T ). 

The least-squares method and the simplex algorithm was 
used (see e.g. Kallrath & Linnell 1987). 

Comparing the theoretical position on the sky Oo and po 
with the observed ones 9i and pi, one can calculate the sum 
of normalized residuals squared 

N 
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Pi - Po,i 



(1) 



following Torres (2004). With this xistr an d using the sim- 
plex algorithm one can obtain a set of parameters (a, p%, i, 
e, 0J3, fl, To) describing the astrometric orbit. 

The weighting is provided by the uncertainties erg and 
a „. These values are obtained from the observations, or es- 
timated as some typical uncertainty level for the certain kind 
of measurement provided by a specific instrument. 

The efficiency and the computing time required by the 
algorithm strongly depends on the initial set of parameters 
and the chosen trial step in the parameter space. If nothing 
is known about the solution, one has to scan a wide range 
of parameters (eccentricity e from to 1, and the angular 
parameters from 0° to 360°). 



The efficiency of the algorithm could be improved if the 
simplex is used repeatedly. It can happen that the simplex 
converges into a local minimum while the global one is far 
away. It is therefore advisable to run the algorithm again, 
with as large initial steps as in the previous run, but keep- 
ing the values of the parameters corresponding to the pre- 
viously found minimum as one vertex. Repeating this strat- 
egy several times over the whole parameter space, one can 
judge whether the global minimum was found by checking 
whether the sum of squares of the residuals is still changing 
or not. 

To conclude, using this strategy and the combined me- 
thod described in section l2~3l one gets the satisfactory result 
after large number of iterations. This number and comput- 
ing time is strongly depended on the input data set and the 
number of parameters fitted. For the case VW Cep with the 
largest data set (about 1600 data points, see below) and also 
with the most parameters to fit (14 in total) one reaches the 
solution, when the sum of squares is not changing signifi- 
cantly, after circa 100,000 simplex steps. This takes about 
one day on computer with 2 GHz processor. 

2.2 Light-time effect 

A different method to study eclipsing binaries in hierar- 
chical triple systems is based on the eclipse timings. The 
light-time effect (or the Tight-travel time') causes apparent 
changes of the observed binary period with a period corre- 
sponding to the orbital period of the third body. This useful 
method is known for decades and its detailed description 
was presented by Irwin in 1959, while some comments on 
it and its limitations were discussed by Frieboes-Conde & 
Herczeg (1973) and by Mayer (1990). 

From the numerical point of view the method is quite 
similar to the previous one, because it is also an inverse 
problem. One has M measurements of the times of min- 
ima of the system at certain constant JDi with the indi- 
vidual uncertainties a m . The task is to find five parameters 
describing the orbit of the third body in the system: the pe- 
riod of the third body p%, the semiamplitude of the LITE 
A, the eccentricity e, the time of the periastron passage T , 
and the longitude of periastron lo 12 . One has to compute 
simultaneously also two (or three) parameters of the eclips- 
ing binary itself, namely its linear (or quadratic) ephemeris 
J Do and period P (and q for the quadratic one). Altogether, 
one has 7 (or 8) parameters to derive from the model fit of 
the minimum-time measurements 

{(JD l ,a m ,i)} i =i,M -> (p 3 ,A,e,T ,uJi2,JDo,P,q). 

Similarly as in the astrometry case, the resultant sum of nor- 
malized square residuals is 

M , {0 _ c) x2 



Xlite 



E 

i=l 



(2) 



http://ad.usno.navy.mil/wds/ 



where O and C stand for observed and computed time of 
minimum, respectively. The same remarks about the effi- 
ciency and the computing time required discussed for the 
astrometric solution also apply here. 
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At this place it is necessary to remark one useful com- 
ment. One has to distinguish between the two angles co^ and 
W12. The parameter used in LITE analysis is W12, but in as- 
trometry the quantity CJ3 = uji 2 + tt is employed. In this 
paper the angle u> stands for the longitude of the periastron 
for the eclipsing binary, i.e. u> — u>i 2 , and the subscripts 
will be omitted for clarity. 

2.3 Combining the methods 

The task is to combine the astrometry and the analysis of 
times of minima into one joint solution together. Having 
7Y astrometric and M minimum time measurements, one 
is able to merge them together and obtain a common set of 
parameters 

(U, 0i, pi,(Te,i , <rp,i, JDi, cr m ,i) — ► 

-> (a,p 3 ,i, e,LO,D,,T , JD ,P,q). 

This set of 10 parameters fully describes the orbit of the 
eclipsing binary around the common centre of mass with the 
third unresolved component together with the ephemeris of 
the binary itself. It is necessary to solve one least-squares fit 
of these 10 parameters. 

One is able to find out the mass of the third body and the 
semimajor axis of the wide orbit because the inclination is 
known and one can calculate the mass function of the wide 
orbit 

(ai2sini) 3 (M 3 sini) 3 
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V 1 — e 2 cos 2 u ' 
where a\ 2 stands for the semimajor axis of the binary orbit 
around the common centre of mass and Mi, M 2 , M3 are the 
masses of the primary, secondary, and tertiary component, 
respectively. For more details see e.g. Mayer (1990). 

The only difficulty which remains unclear is the con- 
nection between the angular semimajor axis a and the am- 
plitude of LITE A. The quantity an could be derived from 
Eq.[3]and with the masses of the individual components one 
is able to calculate also the value 03, i.e. the semimajor axis 
of the third component around the barycentre of the system 
Mi + M 2 



0,3 = a 12 



M 3 



(4) 



The total mutual distance of the components is a to tai = 
Oi2 + a 3- Using the Hipparcos parallax tt (Perryman & ESA 
1997) one can obtain the distance d to the system. Now it 
is possible to enumerate the angular semimajor axis a as a 
function of d and a to tai 

. / Ototal \ 

a = arcsm ^ — - — j . (5) 

The way how these two different approaches were putted 
together is following the similar approach by Torres (2004). 
From the mathematical point of view both methods are anal- 
ogous and there is an overlap of the parameters in both 
methods. Our task is to minimize the combined \ 2 



Xcomfc 



Xlstr 



Xlitei 



(6) 



where the values xistr an d Xlite are taken from the EqsQ] 
and|2] 

Sometimes there arises a problem with the \ 2 values 
in both methods, which are incomparable. Especially when 
there are much more data points in one method against the 
other, also the resultant \ 2 would be much larger and as a 
consequence this method outweigh the other one. This prob- 
lem could be eliminated using new uncertainties a' 

o' B = VN ■ <Jg, <j'=VN-(Tp, <j' m = VIl ■ <J m . 



3 An application of the method to several 
particular systems 

The method described above was tested on a few systems 
satisfying the following conditions: 1. More than 10 times 
of minima and more than 10 astrometric observations are 
available. 2. The observed range of the position angle in the 
astrometric measurements is larger than 10°. 

3.1 VW Cep 

The eclipsing binary VW Cep (HD 197433, BD +75 752, 
HIP 101750, V = 7.3 mag, sp G5V + G8V) is a W UMa 
system and in fact it is one of the most often observed and 
analysed system. Both components are chromosphericaly 
active. VW Cep is rather atypical, because during the pri- 
mary (the deeper one) eclipse is the less massive star (the 
hotter one) occulted by the larger companion (the more mas- 
sive and the cooler one). 

The first observations of its light variations were done 
by Schilt (1926). Since 1946, a large amount of photoelec- 
tric observations was obtained. However, the observed min- 
ima times did not fit the ephemeris due to the LITE and the 
mass transfer between components. There were many light- 
time effect studies of this system and Herczeg & Schmidt 
(1960) proposed the presence of a third body with an or- 
bital period of 29 years and an angular distance of the third 
component between 0.5 " and 1.2 ". 

In 1974, the first successful visual observation of the 
third component was obtained and since then, there were 
16 observations of it. Regrettably, the observations near the 
periastron passage are missing (the gap in data is from 1991 
to 1999). The last two measurements (0 = 231.4°, p = 
0.702" and = 232.6°, p = 0.695") were not published 
yet. These were obtained by a speckle camera in April 2007 
and were kindly sent by Elliot Horch (priv.comm.). 

The most complete set of times of minima is in the most 
recent period study of VW Cep by Pribulla (2000). The first 
times of minima are from the 1920's and altogether 1907 
minima were collected. From this set of times of minima 
313 measurements were neglected due to their large scatter 
(mostly the visual ones). This new minimum-time analy- 
sis is based on a larger data set (about 750 times of min- 
ima more than were used by Pribulla et al.). Two new CCD 
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Table 1 The final results: Parameters of the third bodies from combined astrometry and the light-time effect. The table 
is divided into three parts, in the first one are ten computed parameters, in the second one the values from the literature and 
in the last one the quantities computed from the previous parts. The values of parallax and distance were adopted from the 
Hipparcos measurements. In the row 'Data set', 'a' and 'm' denote the number of astrometric observations, and times of 
minima, respectively. 



Parameter 


Unit 


VWCep 


(Phe 


HT Vir 


J Do 


[HJD] 


2437001.4327 ± 0.0025 


2441643.7382 ± 0.0008 


2452722.5040 ± 0.0050 


P 


[d] 


0.278315349 ± 0.00000012 


1.6697772 ± 0.0000013 


0.4076696 ± 0.0000025 


P3 


[yr] 


29.79 ± 0.08 


220.9 ±3.5 


260.7 ± 0.40 


To 


[HJDl 

L J 


2450390.6 ± 37.3 


2419908.9 ± 2465.6 


2442832.3 ± 60.6 


u> 


[°] 


239.28 ± 2.88 


97.1 ± 2.2 


250.9 ± 0.6 


e 




0.633 ±0.007 


0.366 ±0.082 


0.640 ± 0.005 


A 


[d] 


0.0131 ±0.0004 


0.0808 ± 0.0080 


0.1274 ± 0.0024 


a 


[mas] 


445.2 ± 33.0 


859.5 ± 137.2 


1023.5 ± 183.3 


i 


[°] 


33.6 ± 1.2 


64.4 ± 3.0 


45.4 ±3.7 




[°] 


17.7 ±3.1 


33.5 ±4.9 


180.8 ± 2.6 


M12 


[M ] 


1.37 


6.48 


2.3 


References 




Kaszas et al. 1998 


Andersen 1983 


D'Angelo et al. 2006 


7T 


[mas] 


36.16 ±0.97 


11.66 ±0.77 


15.39 ± 2.72 


d 




27.7 ±0.7 


85.8 ±5.7 


65.0 ± 11.5 


ai2 


[AU] 


4.33 ± 0.20 


15.52 ± 1.60 


32.23 ±0.74 


/(Af 3 ) 


[M ] 


0.0154 ±0.0016 


0.0562 ± 0.0167 


0.1692 ±0.0098 


M 3 


[M ] 


0.74 ± 0.07 


1.73 ±0.26 


2.10 ±0.10 


Data set 




16a + 1594m 


12a + 36m 


275a + 31m 



observations of the minimum light of VW Cep were ob- 
tained in Ondfejov observatory with the 65-cm telescope 
and Apogee AP-7 CCD camera, 2 seconds exposure time 
in R filter. These new times of heliocentric minima are the 
following: 2454154.3094 ± 0.0001 and 2454195.49866 ± 
± 0.00012. 

The short-term variations with the period of about two 
years (see e.g. Kwee 1966 and Hendry & Mochnacki 2000) 
are probably caused by the surface activity cycles on the 
primary component. Due to this activity an unique interpre- 
tation of the behaviour of period changes is still missing. 
Pribulla et al. proposed a mass transfer (the quadratic term) 
plus the third and the fourth body in the system (two peri- 



odic terms). Nevertheless, they were not able to explain the 
O — C diagram in detail. 

Another approach has been chosen in this paper. Espe- 
cially due to only a few astrometric observations (16 mea- 
surements from 1974 to 2007) we have decided to explain 
only the most significant effects in the O — C diagram. It 
means the astrometric variation with a period of about 30 
years has been identified with the O — C variation with 
the same period, but the long-term variation in the O — C 
diagram (which has been interpreted by Pribulla et al. as 
a quadratic term) was explained as a variation due to the 
fourth body on its very long orbit. This approach was chosen 
especially because of the systemic-velocity variations, see 
below. It means during the computation process the value 
Xcomb was minimized with respect to 14 parameters in to- 
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Fig. 1 An O — C diagram of VW Cep. The individual 
observations are shown as dots (primary) and open circles 
(secondary), the small ones for visual and the large ones for 
CCD and photoelectric observations. The curves represent 
the predicted LITE3 + LITE4 (the solid one) and only LITE4 
(the dashed one). 



2425000 2430000 2435000 2440000 
0.06- 



0.04 

5 0.02 
1 

o 

6 -0.02 
-0.04 
-0.06 




-40000 -30000 -20000 -10000 



10000 20000 30000 40000 50000 60000 
Epoch 



Fig. 2 An O — C diagram of VW Cep after subtraction 
of the LITE4. The description is the same as in the previous 
O — C figure, and the solid line represents the LITE caused 
by the third component in the system. 
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Table 2 VW Cep: parameters of the fourth-body orbit 



Parameter 


Unit 


Value ± Error 


Pt 


[yr] 


77.46 ± 0.04 


A 4 


[dj 


0.100 ± 0.002 


U!4 


n 


283.42 ± 2.30 




[HJD] 


2453857.4 ± 13.6 


e4 




0.543 ± 0.007 



„-100 

to 

LU 



tal: 

(A,p 3 , i, e, w, fl, T , JD , P. A i ,p 4 , e 4 , w 4 , T 0)4 ). 

The combined approach of analysing the times of min- 
ima together with astrometry led to the parameters shown 
in Table [TJ and [2] The O — C diagram in Fig. Q] shows the 
times of minima together with the curve which represents 
the LITE3 + LITE4. If one subtracts only the LITE4 varia- 
tion and try to describe the behavior of the minimum times, 
one gets Fig. where only the LITE3 caused by the third 
component is displayed. The fit is not very satisfactory be- 
cause of the presence of the chromospheric activity of the 
individual components, or due to a putative additional com- 
ponent (see e.g. Pribulla 2000). In Fig. [3] the astrometric 
orbit of the binary with the individual measurements and 
their theoretical positions is shown. Regrettably, no obser- 
vations near the periastron passage are available. The curve 
also represents the theoretical orbit according to the param- 
eters given in Table[T]in agreement with the LITE analysis. 

The parameters describing the LITE3 and LITE4 vari- 
ations are in Table Q] and [2] and could be compared to the 
parameters derived during the previous analysis by Pribulla 
et al. Their values for the third-body orbit are: P3 = 31.4 yr, 
e = 0.77, u = 183°, and a tota l = 12.53 AU. Our 
values are in Table [T] except for a to tai = 12.35 AU, and 
as we can see they differ significantly in several parameters. 
This is due to completely different approach describing the 
O — C variations. Only the period and the amplitude of such 
variation are comparable, but these are the most important 
for our combined solution. One has also to disagree with the 
result by Pribulla et al., that the astrometric orbit could not 
be identified with the LITE3 variation from the O — C dia- 
gram. As one can see, our new results are in agreement with 
each other without any problems. 

Also the astrometric orbit could be compared with the 
previously published one. Most recently Docobo & Ling 
(2005) published the following parameters of the astromet- 
ric orbit: p$ = 31.0 yr, a — 485 mas, i = 39.3°, and 
e = 0.68. If one compares these values with the new ones 
(see Table [TJ, one can see that the differences are slightly 
beyond the limits of errors. 

If we assume the mass of the eclipsing binary to be 
M12 = 1.37 Mq (Kaszas et al. 1998) and the parallax 
7r = 36.16 mas (from Perryman & ESA 1997), the dis- 
tance to the system is only about 27.66 pc, which results in 
the third-body mass of M3 = 0.741 M Q . The distant com- 
ponent is about 2.2 magnitudes fainter than the VW Cep it- 
self, so its luminosity and also mass should be much smaller 
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Fig. 3 Relative orbit of VW Cep on a plane of the sky. 
The points represent individual observations, while the solid 
curve corresponds to the solution described in the text. The 
straight lines connect individual observations with their ex- 
pected positions on the fitted orbit. The cross indicates the 
position of the eclipsing binary on the sky and the arrow 
the direction of the orbit of the third body. The dashed-line 
represents the line of apsides. 

than the mass of the eclipsing components. The total bolo- 
metric magnitude of VW Cep is about 4.7 mag, so the mag- 
nitude of the third component is about 6.9 mag, which leads 
to the spectral type of about K3. The typical mass of this 
spectral type is about 0.75 M (according to Harmanec 
1988), which is in good agreement with the new result and 
within its error limits. 

Different systemic velocities v 1 were found at different 
epochs. These values are: u 7 = (—35.4 ± 10) km • s _1 
(Popper 1948), (+9.8 ± 7) km • s^ 1 (Binnendijk 1966), 
(-8 ± 1) km • s" 1 (Hill 1989), and (-16.4 ± 1) km • s" 1 
(Kaszas 1998). In the time plot (see Fig. [4]i one can see 
the curve which represents the theoretical variation of i> 7 



_ 20- 




u -20- 

| -30- 
£-40- 



1930 1940 1950 1960 1970 1980 1990 2000 2010 
Year 

Fig. 4 Systemic velocity variations in VW Cep. The indi- 
vidual points represent computed systemic velocities (see 
details in text). The solid curve stands for the combined 
LITE3 + LITE4 variations, while the dashed one only for 
the LITE4 variation. 
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caused by the orbital motion around the common barycen- 
tre. Except for the first one data point by Popper the sys- 
temic velocities are following the long-term variation and 
are almost within its errors near the theoretical values. The 
value of Popper is affected by a relatively large error. The 
scatter of the individual RV data points is much larger than 
that from Binnendijk, which could be caused by the combi- 
nation of two different data sets from different instruments 
and obtained after more than 600 orbital revolutions (which 
could shift the ephemeris). Pribulla & Rucinski (2006) sug- 
gested that the scatter of the systemic velocity data points 
is instrumental, which seems unlikely for such a large am- 
plitude. For the final confirmation of v 7 variations a more 
accurate and larger data set is necessary. 

We also tried to derive the parallax of VW Cep using 
this combined approach. Leaving the parallax as another 
free parameter, it was calculated from the comparison of the 
angular and absolute semimajor axis. Using this approach, 
almost all of the relevant parameters remained nearly the 
same, only the inclination changed a bit, being about 1° 
lower and the angle £1 about 1.5° lower. The new param- 
eters led to a higher third mass of M3 = 0.76 M©. The 
parallax decreased from (36.16 ± 0.97) mas {Hipparcos) 
to (35.85 ± 0.37) mas. This parallax would shift the dis- 
tance from (27.7 ± 0.7) pc {Hipparcos) to (27.90 ± 0.29) 
pc. As one can see, the result of the parallax is more precise 
than any of the previously derived parallaxes (see Hendry & 
Mochnacki (2000) for a summary of the previous values). 

One could conclude that the third body is probably of 
spectral type K3 with a mass around 0.74 M . It is clear, 
however, that a more complicated model will be needed to 
describe the observed changes completely. 

3.2 C Phe 

The system £ Phe is the brightest eclipsing binary with two 
components of early spectral types, exhibiting total and an- 
nular eclipses. This is the only binary with an eccentric orbit 
included in this study. C Phe (HD 6882, HR 338, HIP 5348, 
V = 4.0 mag, sp B6V + B9V) is an Algol-type eclipsing 
binary. It is a visual triple and SB2 spectroscopic binary. 

The brightest component is the EB, while the most dis- 
tant component is the faintest one (some 6"away and with 
its apparent brightness of about 8 mag). The last known 
component is a 7 th -magnitude star at a distance of about 
600 mas. This is the astrometric component and in our opin- 
ion also the one which causes the LITE. The depth of the 
primary minimum of the eclipsing pair is about 0.5 mag and 
the period is about 1.66 d. 

The unfiltered light curve was observed in 1950's by 
Hogg (1951), after then by Dachs (1971) in UBV filters, 
and the best one by Clausen et al. (1976) in ubvy filters. In 
this latter paper all relevant parameters of the eclipsing sys- 
tem were derived and also the third light was computed. Its 
value changes from 3%(it) to 8%(y) and the distant com- 
ponent was classified as a spectral type A7 star. 
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Fig. 5 The O — C diagram of C, Phe. The description is the 
same as in the previous O — C figures. The apsidal motion 
curve (the dash-dotted one for the primary and the dashed 
for the secondary) is plotted around the (solid) LITE curve. 

Clausen et al. (1976) also collected the times of minima 
obtained before 1975. They concluded that no significant 
apsidal motion is observed. The first apsidal-motion study 
was published by Gimenez et al. (1986). With an updated 
list of the times of minima one is able to conclude that the 
apsidal motion is definitely present. It is clearly visible in 
the O — C diagram shown in Fig. [5] Altogether 36 times of 
minima used here came from the paper cited above and from 
Mallama (1981), Gimenez et al. (1986), Kvfz et al. (1999). 
The most recent ones are in Table [3] 

Our new photoelectric UBV observations were secured 
with the modular photometer utilizing Hamamatsu EA1516 
photomultiplier on the 0.5-m telescope at the Sutherland site 
of the South African Astronomical Observatory (SAAO) 
during two weeks in September 2005. The Johnson UBV 
photoelectric measurements were secured with 10-second 
integration times. Each observation of Q Phe was accompa- 
nied by an observation of the local comparison star -q Phe (V 
= 4.36 mag). All measurements were carefully reduced to 
the Cousins E-region standard system (Menzies et al. 1989) 
and corrected for differential extinction using the reduction 
program HEC 22 rel. 14 (Harmanec & Horn 1998). The 
standard errors of these measurements were about 0.008, 
0.006, and 0.005 magnitude in U, B, and V filters, respec- 
tively. 

The new times of primary and secondary minimum and 
their errors were derived using a least-squares fit to the data 
and by the bisecting-chord method. Only the bottom parts of 
the eclipses were used. The mean values and the errors for 
each individual filter are given. Six new times of minimum 
light were derived using the Hipparcos photometry (Per- 
ryman & ESA 1997) and fitting the published light curve. 
These new times of minima are also included in Table [5] In 
this Table, the epochs are calculated from the light elements 
given in Table Q] the other columns being self-explanatory. 

( Phe has one of the shortest apsidal motions among the 
eclipsing binaries (see e.g. Claret & Gimenez 1993). Due 
to a low eccentricity, the amplitude of the effect is small. 
For an accurate calculation of the apsidal motion rate the 
method described by Gimenez & Garcia-Pelayo (1983) was 
routinely used. The eccentricity of the orbit in the eclips- 
ing binary is e' = 0.0107 ± 0.0020, the longitude of peri- 
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Fig. 6 Relative orbit of £ Phe on a plane of the sky, for a 
detailed description see Fig. [3] Two measurements (the open 
circles) were neglected. 

astron ooq = 12.96° ± 5.96°, and the apsidal motion rate 
Co = (0.028 ± 0.001)°/cycle = (6.16 ± 0.20)°/yr, i.e. 
the apsidal motion period U = 58.5 yr. The most recent 
apsidal-motion analysis is more than 20 years old, made by 
Gimenez et al. (1986), but with no LITE and with a smaller 
set of times of minima. The eccentricity derived by Gimenez 
et al. was almost the same, but the apsidal motion rate Co was 
0.0373 °/cycle and the angle lo = 13°. 

Our approach was a combination of the two different ef- 
fects. The behaviour in the O — C diagram was described as 
a sum of apsidal motion and LITE contribution (O — C) — 
(O - C) apsld + (0- C)li TE , distinguishing the primary 
and secondary minima. It means the least-squares algorithm 
was minimizing the Xcomb w i tn respect to 12 parameters in 
total (A,p 3 , i, e, w, ft, T , JD Q , P, Co, lo q , e'). 

Table 3 The new minima timings of (, Phe based on pho- 
toelectric observations. 



HJD-2400000 Error Prim/Sec Epoch Ref. 



47872.7742 


0.005 


Sec 


3730.5 


[1] 


47873.6172 


0.005 


Prim 


3731.0 


[1] 


48187.535 


0.005 


Prim 


3919.0 


[1] 


48397.0807 


0.005 


Sec 


4044.5 


[1] 


48484.7523 


0.005 


Prim 


4097.0 


[1] 


48508.9557 


0.005 


Sec 


4111.5 


[lj 


51466.9675 


0.0005 


Prim 


5883.0 


[2] 


51467.7907 


0.001 


Sec 


5883.5 


[2] 


53622.6453 


0.0001 


Prim 


7174.0 


[3] 


53623.4710 


0.0001 


Sec 


7174.5 


[3] 



Ref.: [1] - Perryman & ESA 1997; [2] - Shobbrook (2004); [3] - 
This paper; 



The astrometric solution based on the combined app- 
roach is satisfactory, while the older measurements have 
larger scatter than the recent ones (the old ones are visual, 
while the modern ones are speckle-interferometric). Two 
measurements were neglected, because of their large scat- 
ter (see Figj6j. Our solution led to the parameters listed 
in Table Q] which could be compared to previously found 
values. Most recently Ling (2004) reported the parameters: 
p 3 = 210.4 yr, e = 0.348, a = 804 mas, i = 61.9°, 
ft — 33.5°, 003 = 271.7°. It is evident that the new param- 
eters are in very good agreement with these ones. The new 
values imply the mass function of the distant body /(M3) = 
0.056 Mq. With the masses of the primary and secondary 
component of the eclipsing binary Mi = 3.93 Mq and 
M 2 = 2.55 M (Andersen 1983), one could derive the 
mass of the astrometric third body A/3 = 1.73 M Q . This 
value corresponds to a spectral type around A7, which is in 
excellent agreement with the photometric analysis (Clausen 
et al. 1986). 

3.3 HT Vir 

One member of the visual binary STF 1781 is the eclipsing 
binary system HT Vir (ADS 9019, HD 1 19931, HIP 67186, 
BD+05 2794, V = 7.2 mag, sp F8V). HT Vir is a contact 
W UMa system, with a period of about 0.4 d and the depths 
of minima of about 0.4 mag. Both visual components have 
almost equal brightness. The third component of the sys- 
tem is brighter than the eclipsing binary HT Vir during its 
eclipses and fainter than it during its maxima. 

According to Walker & Chambliss (1985) the distant 
astrometric component was discovered by Wilhelm Struve 
in 1830 at a separation of about 1.4 " and position angle 
240°. Since then, numerous astrometric observations were 
obtained (altogether 277, from which 275 were used in our 
analysis) and the orbit is almost completely covered by the 
observations. 

Baize (1972) suggested that the star might be variable. 
After then, Walker & Chambliss (1985) obtained a complete 
light curve of HT Vir and did the first analysis. It indicated 
that both components of the eclipsing pair are almost iden- 
tical and in contact. The temperatures of both components 
are about 6000 K and the spectral type is estimated as F8V 
for the primary and a little bit earlier for the secondary, the 
inclination is close to 90°. The total mass of the eclipsing 
pair is M 12 = 2.3 M Q (D' Angelo et al. 2006). 

Lu et al. (2001) discovered that the distant component is 
also a binary. They have measured the spectra of the HT Vir 
eclipsing pair, and discovered also the lines from the third 
component in the spectra and their RV variations with a pe- 
riod of about 32.45 d. We therefore deal with a quadruple 
system. 

Despite the spectral analysis and a large set of astro- 
metric observations, there were only a few times of min- 
ima published during the last few decades. The main reason 
is relatively recent discovery of the photometric variabil- 
ity of HT Vir. The first times of minima come from 1979. 
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Fig. 7 Relative orbit of HT Vir on a plane of the sky, for 
a detailed description see Fig. [3] 



Since then, there were only 31 observations obtained (see 
Fig. Four new observations were obtained. The two of 
them were carried out in Ondfejov observatory with the 65- 
cm telescope and Apogee AP-7 CCD camera and 1 sec- 
ond exposure time in R filter. This new times of heliocen- 
tric secondary minima are 2454175.58494 ± 0.00011 and 
2454195.56170±0.00007.The next one was observed by L. 
Brat (Private Observatory), using 8-cm telescope with ST-8 
CCD camera, 20 seconds exposure time in R filer, resulting 
in 2454210.44154±0.00015, and the last one by R.Dfeveny 
with ST-7 CCD camera, 60 seconds in R filter, resulting 
in 2454213.49967 ± 0.00012. One unpublished observa- 
tion by M.Zejda was also used and four times of minima 
by M.Zejda published in Zejda (2004) were recalculated, 
because the heliocentric correction was wrongly computed. 

The final plot of the relative astrometric orbit of HT Vir 
is in Fig. [7] The results, parameters of the orbit around the 
common barycentre of the system, are given in TableQ] The 
values of these parameters (A, p$, i, e, u>, £l,To, J Do, P) 
were obtained minimizing the Xcomb- 

Walker & Chambliss (1985) published the first rough 
estimation of the proposed LITE from the astrometric orbit. 
Their value (0.18 d) is not too far from ours (0.13 d). 

The new elements for the astrometric orbit can be com- 
pared to those of Heintz (1986), which are the following: 
p 3 = 274.0 yr, e = 0.638, a = 1010 mas, i = 42.7°, 
fl = 176.4°, uj 3 = 250.0°. As one can see, the period 
of the new orbit is a bit shorter, but the main differences in 
these values are the angles u> and f_. The same fit to the as- 
trometric data points could be reached with simultaneously 
transformed values lo 3 — > C1/3 + 180° and Q — ► f2 + 180°. 
This only means the interchange of the role of the two com- 
ponents. This result indicates the incorrect identification of 
the variable HT Vir in the system in our analysis (the vari- 
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Fig. 8 An O — C diagram of HT Vir. The description is the 
same as in the previous O — C figures, all minimum times 
are photoelectric or CCD ones. 



able was supposed to be the component A) and also in the 
WDS catalogue, see WDS notesB While Pribulla & Rucin- 
ski (2006) correctly identified the variable HT Vir as the B 
component and A as a single-lined binary. 

If we adopt these parameters to estimate the mass func- 
tion of the distant pair (mass function of the whole pair, not 
the individual components), we obtain /(M3) = 0.17 M Q . 
This is quite a high value, dictated by the large amplitude of 
the LITE. With the total mass of the primary and secondary 
M12 = 2.3 M we get the third mass of M 3 = 2.10 M . 
The mass of the distant pair is quite high (D'Angelo et al. 
(2006) derived the mass for some 50 per cent lower, M3 = 
1.15 Mq), but note that also this object is a binary and we 
do not know the individual masses. From the spectroscopic 
observations (we remind that it is a SB1 type binary), we 
are only able to estimate the mass function of the compo- 
nents, or some upper limit for one of them (we do not know 
the inclination). Our resulting M3 is the total mass of the 
SB1 pair M31 + M3 2; the limit for the invisible-component 
mass M3.2 • sin(i) = 0.075 M . If we assume the copla- 
nar orbit, high difference in masses would arise, one com- 
ponent should be much more luminous and also more lumi- 
nous than the eclipsing pair itself, which is not the case. In 
fact the whole system is not coplanar (see e.g. i — 315.5° 
and the inclination of the EB close to 90°). If we assume 
two approximately equal masses, there is a problem with 
the luminosity, because the distant pair has to be roughly as 
luminous as the eclipsing pair. This could only be satisfied 
if one component is underluminous or degenerate. 

We have to take into consideration also the comment on 
the light-curve solution by Walker & Chambliss (1985). Us- 
ing the Wood's model, they discovered that if the third light 
L3 is fixed to be equal to the light from the distant visual 
component (it means L3 = 0.5), the solution of the light 
curve is unrealistic. To conclude, the system could be much 
more complicated than the approach we have used here. Es- 
pecially because of the resultant mass and luminosity of the 
distant pair, the body causing the astrometric variation is 
probably different from the one causing LITE, but this con- 
clusion will be proven only if also the nonlinear part of the 
O — C diagram is covered. 



http://ad.usno.navy.mil/wds/wdsnewnotes_main.txt 
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4 Discussions and conclusion 

Although the number of systems where the astrometric orbit 
of the third component has been measured together with the 
presence of LITE is growing steadily, in the most cases only 
very limited coverage, both in astrometry and times of min- 
ima is available. Especially due to this reason the combined 
analysis of these systems is still difficult, the results are not 
very convincing and the resultant parameters are affected by 
relatively large errors. 

During the last decade a few papers combining the ap- 
proach of simultaneous solution of radial velocities, spec- 
tral analysis, astrometry, Hipparcos measurements or LITE 
were published. Besides the systems mentioned in the in- 
troduction (V505 Sgr, QS Aql, 44 Boo, QZ Car, SZ Cam, 
GT Mus, and V2388 Oph) there were also the analyses of 
VI 061 Cyg (combining the light curve analysis, radial ve- 
locity analysis, light-time effect and Hipparcos measure- 
ments, see Torres et al. 2006), the papers where the radial 
velocity measurements and astrometry were combined (see 
Muterspaugh et al. (2006) for the solution of LITE system 
V819 Her, or Gudehus (2001) for ^ Cas), the paper on HIP 
50796 combining the radial velocities with the Hipparcos 
abscissa data (Torres 2006), or the paper on 8 Lib compar- 
ing the results from the period analysis, light-curve analy- 
sis, spectral analysis, radio emission and astrometry, respec- 
tively; see Budding et al. (2005). 

Three eclipsing binaries were studied in this paper. In 
the case of VW Cep, where the orbits in both methods have 
relatively best coverage, the distant body satisfies the limits 
for the luminosities, and also the systemic velocity varia- 
tions coincide with our hypothesis. New results are compa- 
rable with the previous ones. An additional fourth body was 
introduced to describe the long-term variation in times of 
minima. The system is probably more complicated than we 
assumed (chromospheric activity cycles, stellar spots and 
flares), and we decided to explain only the most pronounced 
effects in the O — C diagram. Using this combined approach 
it is possible to derive the parallax of VW Cep more pre- 
cisely than in previous papers, resulting in it = (35.85 ± 
± 0.37) mas. The system ( Phe displays an apsidal mo- 
tion together with LITE and this explanation fits the O — C 
residuals quite well. The quality of the astrometry is worse, 
but the main result about the masses and the spectral type 
of the third body is in an excellent agreement with the pre- 
vious photometric analysis. The last system HT Vir is the 
case where the new resultant parameters of the distant body 
are about 2 times larger than one would expect. This could 
be due to the fact that only a few times of minima were ob- 
served in the linear part of the O — C diagram, new minima 
are needed in the next decades. 

The final result is that the method itself is potentially 
very powerful but it is also very sensitive to the quality of 
the input data, especially if the method is used for determin- 
ing the distances of these binaries. It can only be applied 
successfully in those cases where the astrometric orbit and 



the LITE in the O — C diagram are well defined by existing 
observations. 
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